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ABSTRACT 

We present line and continuum measurements for 9818 SDSS type-I active galactic nuclei (AGNs) with 
Z < 0.75. The data are used to study the four dimensional space of black hole mass, normalized accretion rate 
(L/LEdd), metalicity and redshift. The main results are: 1. LfL^dd is smaller for larger mass black holes at all 



redshifts. 2. For a given black hole mass, L/L 



Edd c 



or(l 



where the slope y increases with black hole 



mass. The mean slope is similar to the star formation rate slope over the same redshift interval. 3. The FeII/H/3 
line ratio is significantly correlated with L/Lgdd- It also shows a weaker negative dependence on redshift. 
Combined with the known dependence of metalicity on accretion rate, we suggest that the FeII/Hj3 line ratio is 
a metalicity indicator. 4. Given the measured accretion rates, the growth times of most AGNs exceed the age 
of the universe. This suggests past episodes of faster growth for all those sources. Combined with the FeII/H/3 
result, we conclude that the broad emission lines metalicity goes through cycles and is not a monotonously 
decreasing function of redshift. 5. FWHM([Olll] A5007) is a poor proxy of a* especially for high L/LEdd- 
6. We define a group of narrow line type-I AGNs (NLAGNls) by their luminosity (or mass) dependent H/3 
line width. Such objects have L/Lsdd> 0.25 and they comprise 8% of the type-I population. Other interesting 
results include negative Baldwin relationships for EW(Hj3) and EW(Fell) and a relative increase of the red part 
of the Hj3 line with luminosity. 

Subject headings: Galaxies: Abundances - Galaxies: Active - Galaxies: Nuclei - Galaxies: Starburst - 
Galaxies: Quasars: Emission Lines 



1. INTRODUCTION 

Recent progress in reverberation mapping of type-I (broad 
emission lines) active galactic nuclei (AGNs) provide a re- 
liable method for deriving the size of the broad line region 
(BLR) as a function of the optical continuum luminosity, the 
Hj3 line luminosity, the ultraviolet continuum luminosity and 
the X-ray (2-10 keV) luminosity (Kaspi et al. 2005, hereafter 
K05 and references therein). This has been used in numerous 
papers to obtain a "single epoch" estimate of the black hole 
(BH) mass by combining the derived BLR size with a mea- 
sure of the gas velocity obtained from emission line widths, 
e.g. FWHM(HJ3) (e.g. Netzer 2003; Vestergaard 2004; Shem- 
mer et al. 2004, hereafter S04; Baskin and Laor 2005). There 
are obvious limitations to such methods due to the scatter in 
the BLR size, the source luminosity and line width. These 
translate to a factor of 2-3 uncertainty on the derived masses 
when using the H/3 line width as the velocity measure and 
the optical continuum luminosity (XL^ at 5100A, hereafter 
L5100) f° r estimating the BLR size. The uncertainty is larger, 
and perhaps systematic, when other emission lines are used. 
Additional uncertainties are associated with the currently lim- 
ited luminosity range, 10 42 <L 5100 < 10 46 erg s _1 , thus mass 
determination for very luminous sources must be based on ex- 
trapolation. Notwithstanding the uncertainties, such methods 
can provide BH mass and accretion rate in large AGN samples 
that cannot be obtained in any other way. 

This paper follows several earlier studies of type-I AGNs in 
the Sloan Digital Sky Survey (SDSS; see York et al. 2000) 
spectroscopic archive. The main aims are to obtain BH 
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masses and accretion rates for a large number of such sources 
and to understand the time evolution and the metal content 
of the BLR gas. Earlier studies of this type include, among 
others, McLure and Jarvis (2002) and McLure and Dunlope 
(2004). Here we use the latest SDSS public data archive and a 
much larger number of sources compared with earlier studies. 
This enables us to obtain statistically significant correlations 
for various sub-samples of the data. We also introduce new, 
more robust methods for measuring the emission lines thus 
enabling more reliable study of metalicity and accretion rate. 

This paper is organized as follows. In §2 we introduce our 
sample and explain the method of analysis and the various 
measurements performed on the SDSS spectra. In §3 we show 
the more important correlations pertaining to type-I sources 
and in §4 we give a discussion of the main results, including 
the derived redshift dependent accretion rate which we com- 
pare with the known star formation rate over the same redshift 
interval. 

2. BASIC MEASUREMENTS 

2.1. Sample selection 

The present sample contains all broad line AGNs in the 
SDSS sample from data release 4 (DR4, see Adelman- 
McCarthy et al. 2006 and references therein) with a redshift 
cut at z = 0.75. This more than doubles the number of such 
sources in DR3 (Schneieder et al. 2005). The sources are 
defined as QSOs by the target selection algorithm but will be 
addressed here as type-I AGNs. The limit of 0.75 is dictated 
by the wavelength band of the SDSS spectroscopy and is de- 
signed to allow the measurement of the H/3 and [O III] A5007 
lines and the 5100A continuum in all sources. The inclusion 
of [Olll]A5007 is required for different line velocities and 
for omitting very few sources with such weak broad lines that 
they are confused with type-II AGNs. The requirement on 
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H/3 is motivated by the work of Kaspi et al. (2000) and K05 
which present the only available reverberation mapping sam- 
ple. This low redshift sample includes H/3 and L5100 measure- 
ments for 35 sources and much fewer UV line and continuum 
measurements. The BLR size is obtained from the measured 
H/3 lag (or a combination of Ha and H/3 lags) and the veloc- 
ity from the observed FWHM(Hj3). Different lines have been 
used in several other papers (e.g. Vestergaard 2004; Vester- 
gaard and Peterson 2006) yet their use is more problematic 
since it requires a secondary calibration against Hj3 and Lsioo- 
Of the two UV lines most commonly used, Mg II A2798 is 
probably more reliable (e.g. McLure and Jarvis 2002; McLure 
and Dunlop 2004) and C IV A 1549 suffers from various sys- 
tematic effects (Baskin and Laor 2005 and references therein). 
We limit the present study to the use of H/3 and defer the study 
of higher redshift AGNs, based on the Mg II A2798 line, to a 
later publication. 

The SDSS archive (Stoughton et al. 2002) was queried for 
QSOs with z < 0.75. This returned 14554 objects out of which 
1212 were found to have problematic mask data (spectra with 
more than 30% flagged pixels; see Stoughton et al. table 1 1 
for details). These objects are not included in the analysis. 
We decided to stick with the original magnitude limit of the 
sample at the time of target selection. This means that QSOs 
which are fainter than i(mag)=\9.l in the TARGET database, 
but were later added to the BEST database due to improved 
photometry, or objects that were added to the QSO database 
after their spectra were examined as different science targets 
(thus listed as TARGET_QSO_FAINT by the SDSS target se- 
lection algorithm; see §4.8.4 and table 27 in Stoughton et 
al.) were not included. This removes 1927 additional sources 
from the sample leaving 1 1415 sources with i(mag) < 19.1. 

The spectra of all sources were corrected for Galactic ex- 
tinction using the dust maps of Schlegel et al. (1998) and 
the Cardelli et al. (1989) model. Next, pixels with problem- 
atic MASK values (as defined above) are replaced by local 
interpolation and the spectra are de-redshifted and put on a 
uniform wavelength scale of 1 A per pixel. 

2.2. Line and continuum measurements 

Our work make use of the SDSS spectra as obtained from 
the archive without any attempt to subtract the galaxy con- 
tribution. This introduces some uncertainties at the low lumi- 
nosity end of the sample and is addressed below in those cases 
where we consider it to be important for the derived correla- 
tions. The general issue of stellar contribution to AGN spectra 
of various luminosities is discussed in Bentz et al. (2006). 

The emission lines measured in this work are [O III] A5007, 
Hj3 and the Fell line complex between the wavelengths of 
4000 and 5300A. We first fit a linear continuum between two 
continuum bands: 4430-4440A and 5080-5 120A. This con- 
tinuum is subtracted from the spectrum and a five Gaussian 
model is fitted to the profiles of the two [O III] lines and H/3 . 
The [O m]AA4959, 5007 A lines are forced to have their theo- 
retical intensity ratio (3.0), the same FWHM and their known 
wavelengths. The H/3 line is fitted with two broad Gaussians 
and a single narrow Gaussian with FWHM and wavelength 
that agree with the [O III] lines to within 200 km s _1 . At this 
stage, the narrow Hj3 component is not allowed to exceed 0.7 
the intensity of [O III] A5007. The broad H/3 Gaussians are 
limited to FWHM of 1500-5000 km s" 1 and a shift of up to 
500 km s _1 relative to [Om] A5007. These first limits are 
necessary in order to avoid the inclusion of some neighboring 



Fell lines, mostly on the blue side of the Hj3 profile, and to 
properly model the broad H/3 components. 

The next step involves fitting and deblending of Fell lines. 
The Boroson & Green (1992; hereafter BG92) Fe II emis- 
sion template was convolved with a single Gaussian, which 
has the same FWHM as the broad Hj3 from the previous 
fit. The resulting Fell model is fitted to the spectra over 
the region 4400-4650A. This region was chosen to mini- 
mize the effect of contamination with H/3, HeIlA4686 and 
[O III] A5007. Next we measure the contribution of the Fell 
lines to the previously defined continuum bands, subtract this 
and refit the continuum. The entire procedure is repeated 
and the newly defined Fell template is subtracted from the 
continuum-subtracted spectrum. This leaves us with a Fell- 
free spectrum. 

The final stage involves, again, a five Gaussian fit to the 
[O III] and H/3 lines. At this stage we relax the conditions 
and allow a broader range of Gaussians (as narrow as the 
[Om]A5007 line and as broad as 20,000 kms" 1 ), and a 
larger shift between the narrow H/3 and [O III] A5007 (450 
km s~'). The 20,000 km s _1 limit is required since almost 
all cases exceeding this limit are found to show poor quality 
fits (see below). We also allow shifts of up to ±1000 km s -1 
between each of the two broad H/3 components and the wave- 
length of the narrow H/3 line. All line and continuum mea- 
surements used later in the statistical analysis are based on 
the final stage of the fit and the second Fell fit. We found this 
multi-stage method to be more robust, and less prone to fit- 
ting inaccuracies compared with a single-stage fit of all line 
and components. 

Regarding FWHM(fO III] A5007), we use an "intrinsic" 
line width which is obtained by properly subtracting a uni- 
form instrumental line width of FWHM=140 kms~'. We 
have decided not to include in the fit the blue wing commonly 
observed in the line profile. This removes a potential compli- 
cation due to blending with the broad H/3 profile and is also 
in line with earlier measures of the line in type-I and type-II 
sources. As explained in Greene & Ho (2005), a single com- 
ponent [O III] A5007 line is a good indicator of a*. This point 
is further discussed in §3.5. We do not include in the analy- 
sis the He II A 4686 line which is clearly seen in many of the 
sources. 

As for the broad H/3 line, we measure its rest equivalent 
width (EW(Hj3)) and FWHM(Hj3) and also divide the com- 
bined line profile into two parts, at the rest wavelength of 
4862A (as determined by the wavelength of [Om] A5007). 
The line flux is then measured in each part. We note that 
some earlier works followed BG92 and use EW(Fell) as a 
replacement for the iron line intensities. Here we define 
L(FeII) by integrating the line flux over the wavelength range 
of 4434-4684A (e.g. S04). In what follows FeII/H/3 stands 
forL(FeII)/L(H/3). 

The automated procedure to measure line and continuum 
fluxes suffers from various problems and limitations, some 
that are difficult to cure. To allow better control of the fit 
results, we make use of the R 2 and % 2 statistics (both per de- 
gree of freedom). These are combined in various ways to re- 
ject sources that do not fall within the predefined limits of 
R 2 > 0.2 and % 2 < 5. A common problem is encountered in 
weak flux sources with very large FWHM(H/3). In many of 
these cases, mostly those at small redshift, the stellar contri- 
bution to L5100 is of the same order as the non-stellar contin- 
uum. The poor signal-to-noise (S/N) in some of those makes 
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FIG. 1. — Basic sample properties. Top panel: L5100 vs. redshift for RL 
(magenta points) and RQ (black points) type-I AGNs. Middle panel: M(BH) 
vs. redshift for type-I RQ-AGNs. Bottom panel: L/L Edd vs. M(BH) for 
type-I RQ-AGNs. Here and in all other diagrams, small black points mark 
the entire RQ sample while the colored dots represent sources in two mass 
groups: II) 5 s M e (red) and I If' v> M (blue). 



it difficult to properly fit the low contrast Hj3 profile. A sim- 
ilar problem occurs with very broad low contrast Fell blends 
which can result in erroneous, very large FeII/Hj3. Such ob- 
jects are usually found and rejected by their small value of 
R 2 . Given the strong dependence of luminosity on redshift, it 
is not surprising that such problems are less common at higher 
redshifts. Thus, the fit quality at high redshift is usually better 
even in sources that have low S/N due to being close to the 
flux limit of the sample. Our rather conservative approach re- 
sults in the rejection of about 13% of all sources due to poor 
fits. Relaxing those criteria to include about half of the re- 
jected sources in the analysis change very little the conclu- 
sions of the paper. 

Radio information is readily available from the SDSS 
archive that provides FIRST (First Images of the Radio Sky 
at Twenty centimeters) data for all sources (e.g. Richards et 
al. 2002; McLure and Jarvis 2004). This enables us to divide 
the sample into radio loud (RL) and radio quiet (RQ) AGNs. 
The final sample that it analyzed here includes 9818 sources 
out of which 812 are RL-AGNs. Its luminosity distribution as 
a function of redshift is shown in Fig^ 

2.3. Luminosity black hole mass and accretion rate 

The combination of L5100 and FWHM(H/3) was converted 
to BH mass using the most recent reverberation mapping re- 
sults of K05. This gives 



M(BH) = 1.05 x10 s 



^5100 



10 46 



erg sec 



0.65 



FWHM(Hp ) 
1000 km /sec 



Mr. 



(1) 

The slope chosen here (0.65) is obtained from the RBLR-L5100 
relationship of K05. More recent papers (Vestergaard and Pe- 



terson 2006; Bentz et al. 2006) suggest a flatter slope (0.52) 
based on improved host galaxy subtraction in several of the 
low luminosity sources in the K05 sample, and the omission 
of several sources from the original list. We note that the 
inclusion of all the K05 sources in the Bentz et al. (2006) 
analysis, using the BCES method (see the Bentz et al. paper 
for explanation and their Table 1 for results) gives a slope of 
0.57 ±0.08 which is entirely consistent with the K05 slope of 
0.66 ± 0.07. Moreover, there are other reasons for choosing a 
slope steeper than 0.52, e.g. using L(Hj3) as a proxy for the 
UV continuum luminosity (which is not affected by the host 
galaxy subtraction) gives such a slope. Finally, removing sev- 
eral of the faintest sources from the K05 sample also justifies 
this slope. For example, in 91% of the sources in our sam- 
ple L5ioo> 10 435 erg s -1 . Thus it is justified not to include 
the K05 sources with Lsioo< 10 435 erg s _1 when searching 
for the most suitable reverberation mapping relationship for 
the SDSS sample (this also avoids complications due to host 
galaxy contamination). This smaller sample produces a rela- 
tionship very similar to the one given by eqn.^ 
The normalized accretion rate, L/Lem, is given by 

r it _ /zXsioo 

1 Edd ~ 1.5 x 1O 38 [M(BH)/M ] ' K ) 

where fi is the bolometric correction factor. Estimates of /l 
requires further assumptions about the spectral energy distri- 
bution (SED), in particular the luminosity of the unobserved, 
extreme UV continuum. This issue has been discussed in 
great detail in numerous publication. Estimates for fi in type- 
I AGN range between 5 and about 13 (e.g. Elvis et al. 1994; 
Kaspi et al. 2000; Netzer 2003; Marconi et al. 2004; S04). 
Excluding dust reprocessed radiation (that was included in the 
Elvis et al. 1994 SED), the range is about 5-10. There are 
clear indications from general considerations of the AGN lu- 
minosity function (Marconi et al. 2004), as well as from direct 
spectroscopy (e.g. Scott et al. 2004; Shang et al. 2005), that 
the shape of the SED is luminosity dependent and the ratio is 
larger in lower luminosity AGNs. Since the exact luminosity 
dependence is not known (see however Marconi et al. 2004, 
Fig. 3 for a specific suggestion) we adopt a constant ratio of 
fi = 7 which is roughly the mean between the value deduced 
for weak and strong UV bump sources. A more advance anal- 
ysis can involve the actual Marconi et al. values. For the range 
ofL, 

ioo considered here (Fig. 1), this would cover the range 
6.2 < / L <8.3. 

The rest of this paper includes various diagrams involving 
mass, luminosity and accretion rate of AGNs at different red- 
shifts. Despite the specific expressions given in eqns.^and 
|2j none of those correlations is trivial since the dependence 
of luminosity on the BH mass is unknown a priori. Our work 
can be considered as a study of the four dimensional distribu- 
tion of type-I AGNs along the coordinates of M(BH), M(BH), 
metalicity and redshift. 

The SDSS is a flux limited sample and hence many of the 
derived correlations are affected by its redshift dependence 
incompleteness. In the following we take into account the 
flux limit of the sample which is transformed, as appropri- 
ate, to the various correlations under study. The limits on 
L5100, L/LEdd etc. are obtained from the (reddening cor- 
rected) i(mag), the redshifts and a k-correction that assumes 



-0.5 



across the wavelength range of interest (5100- 



8000A). 

Fig[^shows BH mass as a function of redshift for all RQ- 
AGNs in our sample. It confirms the well known tendency for 
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finding larger active black holes at higher redshifts, a result 
that was already presented by McLure & Dunlope (2004) over 
a wider redshift range. The lower boundary at each redshift 
is a clear manifestation of the strong correlation of M(BH) 
and L5100 in our flux limited sample. The diagram is very 
similar to the one presented in McLure & Dunlope (2004) 
where a smaller sample and somewhat different expressions 
for deriving M(BH) and L/LEdd have been used. Here, and 
in several other diagrams, we focus on two mass groups that 
are marked in two different colors: 10 7 ■ 5 ~ 7 - 8 M© (red, 1207 
sources) and l0 8 - 5 ~ 8 - 8 M (blue, 1422 sources). 

McLure and Jarvis (2004) used a sample of about 6,000 
SDSS QSOs to investigated the mass of the central BHs as 
a function of radio properties. They found that RL-QSOs 
host, on the average, more massive BH for a given optical 
continuum luminosity. The mean M(BH) difference between 
RL-AGNs and RQ-AGNs found by those authors is 0. 16 dex. 
This was attributed to the larger mean FWHM(Hj3) in radio 
loud sources. We checked these findings by comparing the 
luminosity, BH mass and L/LEdd for the two populations in 
our z < 0.75 sample. Using the medians of all properties, 
we find a difference of 0.11 dex in L5100, a factor of 1.06 in 
FWHM(HJ3), and a difference of 0.15 dex in M(BH) (the dif- 
ferences in the mean are somewhat larger). The difference in 
M(BH) is, of course, the result of the differences in luminosity 
and line width between the two populations (see eq. Q). 

Fig.[T]also shows L/L Edd vs. M(BH) for the RQ-AGN sam- 
ple as well as for the two mass groups. The low accretion rate 
limit is determined, mostly, by the flux limit of the sample. 
Thus, we expect that more object with small L/L^dd are miss- 
ing at higher redshifts. The highest accretion rate sources in 
each mass group are not affected by this limit. Based of these 
sources we confirms that the lowest mass BHs are the fastest 
accretors, a result that was already noted by McLure and Dun- 
lop(2004). 

3. RESULTS 

3.1. Accretion rate vs. redshift 

The correlation of L/LEdd with redshift for the entire sam- 
ple, as well as for the two mass groups, is shown Fig.|2] The 
diagram shows also two computed lower limits for each of the 
two groups corresponding to the flux limits of the two mass 
boundaries. The visual impression is somewhat misleading 
since much of the apparent correlation is due to the lower limit 
on L/LEdd resulting from the flux limits on the various mass 
groups. The presence of such flux limits considerably com- 
plicate the analysis and we have searched for rigorous ways 
to derived the intrinsic z-L/LEdd dependence. Below we de- 
tail the two methods that were used, the peak of the L/L Ed d 
distribution and the Efron and Petrosian (1992) truncated per- 
mutation test. 

3.1.1. The peak distribution method 

We first study the distribution of L/L Ed d m several redshifts 
bins that are just wide enough to include a reasonable number 
of sources. These histograms (see Fig. [3} are used to define 
the peak and the envelope of the L/LEdd vs - redshift relation- 
ship. The bins are defined such that Az/z = 0.1 and include 
from 50 to 230 sources, with an average number of about 120 
sources per bin. 

Such a procedure is only meaningful if the peak of the dis- 
tribution is clearly defined, i.e. if the value of L/LEdd corre- 
sponding to the maximum number of sources is larger than 




FIG. 2.— L/L Edd vs. redshift for our SDSS sample of type-I RQ-AGNs. 
The solid lines represent the limiting accretion rates obtained from the flux 
limits for the two mass groups (see text). Symbols are the same as those of 
Fig.l. 
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FIG. 3. — L/LEdd histograms for various redshift bands, as marked, for the 
low (bottom) and high (top) mass groups. The redshift bins are defined in the 
text. The vertical lines at the top of each panel mark the conservative lower 
limit on L/L^dd obtained from the flux limit of the same color histogram (see 
text for more details). 



the accretion rate imposed by the flux limits of the sample. 
There are four such limits for each bin: two corresponding to 
the redshift boundaries and two to the mass boundaries. The 
largest of those four (i.e. the most conservative limits) are 
shown in Fig.[3]for each of the two mass groups and several 
redshift bins. As evident from the diagram, the peaks of the 
L/LEdd distribution in the lower mass group (bottom part of 
the diagram) are only meaningful for z < 0.3. Beyond this 
redshift the method is not applicable since the SDSS survey 
is not deep enough to properly sample many sources and the 
limiting L/LEdd is larger than the value at the peak. The larger 
mass group (top panel) is less problematic and the histograms 
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FIG. 4. — L/LEdd vs. z for various mass groups. The color symbols cor- 
respond to the two mass groups where solid line represent peaks and dotted 
lines represent the upper envelopes of the distributions, as explained in the 
text. Black lines and symbols represent a third mass group with M=10 s8 ~ 7 ' 2 
Mq . The straight solid lines are the results of the truncated permutation test 
for the same mass groups (eqn|4]and Table 1). 



show that the peaks are well defined for 0.3 < z < 0.75. The 
z=0.2 histogram (not shown here) can also be used to define 
the peak. 

As an independent check on the slope of the correlation 
in those redshift intervals where the peak distribution can be 
identified, we also used the upper envelope of the distribution 
which we define by the fifth largest L/L E dd (the number of 
sources with larger L/LEdd ls so small in some groups that 
we suspect incomplete sampling). Obviously, the upper en- 
velop is not affected by the flux limit of the sample. How- 
ever, the redshift bins must contain enough sources to justify 
this method. This may be a problem for very high luminosity 
sources at low redshift because of the limited area covered by 
the SDSS sample. 

The results of the test are shown in Fig. |4] We also show 
the peaks of the distribution for a third mass group, with 
M= l0 6 - 8 ~ 7 2 M , for those redshifts where the peak of the 
L/LecUI distribution is higher than the flux limit of this group. 
In this mass range, this is only true for z < 0.25 since there 
are very few such low mass sources at higher redshifts be- 
cause of the flux limit. The results can be presented as sim- 
ple power-laws in redshift ranging from L/LEdd'* z° ' 8 for for 
M(BH)=1O 7 M to about L/L Ed d^z 22 for M(BH)=10 86 M Q . 
Given the small number of points in the various redshift inter- 
vals, and the differences between the slopes obtained from 
the peaks and the envelopes of the various histograms, we es- 
timate an uncertainty of about 0.3 on the slopes obtained with 
this method. 

3. 1 .2. The truncated permutation test 

The determination of bivariate distributions from truncated 
astronomical data is a well known problem. Such data sets 
are very common and the flux limited sample considered here 
is only one such example. The problem is to test such trun- 
cated bivariate data sets for statistical independence without a 
priori knowledge of the real underlying correlation. This is- 
sue is discussed in great detail by Efron & Petrosian (1992; 
see this paper for earlier attempts to address this issue) who 
suggest powerful nonparametric permutation tests that can be 
applied to data sets of this type. The paper describes tests for 
nontruncated as well as truncated data sets and apply them to 



TABLE 1 

Fit coefficients forL/L £iW vs. z (eqns. 4 & 5) 



logM (M s ) 


Az 


7 (90% limits) 


S (90% limits) 


7.0 


0.05-0.3 


1.02 (0.91-1.13) 


8.6 (7.7-9.7) 


7.5 


0.05-0.75 


1.25 (1.15-1.34) 


7.6(7.1-8.05) 


8.0 


0.1-0.75 


1.61 (1.54-1.67) 


6.8 (6.6-7.1) 


8.5 


0.20-0.75 


1.79(1.70-1.87) 


6.1 (5.7-6.3) 


9.0 


0.25-0.75 2.30(2.15-2.47) 


7.6(7.1-8.1)) 



two astronomical cases that are rather similar to the sample 
in question. In particular, they describe a rank-based statistics 
that is powerful, robust and easily administrated to data sets 
like ours. 

The normalized rank statistic of Efron and Petrosian (1992), 
Tj, is used to define a test statistic 

f lv ( fl f fl f fl )=E i w i ^/(E i w 2 ) 1 / 2 , (3) 

where w= (w\ ,W2, ■■■■w„) is a vector of weights. t w (data) has 
mean and variance 1 and is distributed normally thus allow- 
ing to obtain the best value and the confidence limits using 
the normal approximation (see Efron & Petrosian 1992 §3 for 
more details and definitions). The choice made here, and rec- 
ommended also in Efron & Petrosian (1992), assumes equal 
weights for all points, i.e. w= (1,1, ....1). The application of 
the permutation test involves a guess of the real distribution 
which recovers the null hypothesis of independence between 
the two variables in questions (in our case L/LEdd an d z )- 

We have used truncated permutation tests to check two fam- 
ilies of ad hoc luminosity redshift dependences: 

L/LEdd -z y(M) , (4) 

and 

L/L E dd~{l+z) S(M) . (5) 

This involved changing the guessed values of y and 8, for var- 
ious mass groups, and looking for their best values (i.e. those 
7 and 8 resulting in t w — 0) as well as the 90% confidence 
limits (those value giving t w = —1.645 and t w = 1.645). The 
mass groups used were limited to be 0.2 dex wide except for 
the lowest mass group, around M(BH)=10 7 M , where the 
small number of objects dictates a wider range (0.4 dex). The 
results are listed in Table 1 where we also specify the redshift 
range for each group (the largest mass group contains no ob- 
jects at small z and the lowest mass group contains very few 
sources at z > 0.25 hence the need to define a redshift range). 

The results of the test are shown in Fig. |4] for the case of 
eqn. |4] They are consistent with the results of the peak dis- 
tribution method and lend support to the interpretation that 
the peak of the distribution has indeed been recovered for the 
mass groups and redshift intervals discussed above. These re- 
sults are further discussed and compared with the cosmic star 
formation rate in §4.2 

3.2. Fell/Hp and iron abundance 

Fig-IUshows the dependence of FeII/Hj3 on L/LEdd for our 
sample of type-I SDSS AGNs. As clear from the diagram, and 
has been verified by standard regression analysis, FeII/H/3 
and L/LEdd are significantly correlated and larger L/LEdd are 
associated with stronger Fell and/or weaker Hj3 lines. This 
conclusion applies for all sources with L/LEdd ~ 0.03 while 
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FIG. 5.— FeII/H/3 vs. L/L £(W for our SDSS sample of type-I RQ-AGNs. 
Symbols are the same as those of Fig. 1 . 



at lower accretion rates the scatter is too large to see a definite 
trend. Such a correlation has been suspected in the past and 
has studied in several previous papers (see Netzer et al. 2004 
and references therein). However, it was never tested in such 
a large sample and over the entire 0-0.75 redshift range. 

We have also studied the dependence of FeII/H/3 on red- 
shift in the entire sample as well as in the various mass 
groups. Simple linear regression suggests that, for RQ-AGNs, 
FeII/Hj3oc z~° ' 22±0 - 02 . The slope is somewhat steeper for RL 
sources (—0.38 ±0.07) and for the mass group of l0 7 - 7 ~ 7 - 8 
Mq objects (—0.40 ± 0.05). The large number of points in all 
those groups of sources results in formally significant corre- 
lations. However, detailed examination shows that the results 
depend strongly on the measured FeII/H/3 at the low redshift 
low luminosity end of the distribution. As explained in §2.2, 
some of these cases, especially those with extremely broad 
lines, are difficult to measure and are therefore uncertain be- 
cause of the large stellar contribution. There are indeed dif- 
ferences in the mean FeII/H/3 between the z < 0.2 and z > 0.5 
redshift bins but not a smooth continuous change as a func- 
tion of redshift. In view of this, we consider the evidence for 
redshift dependence of FeII/H/3 as marginal. This uncertainty 
does not affect the FeII/H/3 vs. L/LEdd correlation discussed 
earlier. 

The modeling of Fell line formation in AGN clouds is ex- 
tremely complicated. It involves hundreds of energy levels 
and thousands of transitions. Earlier detailed studies by Net- 
zer and Wills (1983) and Wills, Netzer & Wills (1985), that 
included most of the important atomic processes, reached the 
conclusion that the Fell line intensities are too uncertain to 
be used as iron abundance indicators. Major improvements 
in atomic data (Sigut & Pradhan 2003), combined with im- 
proved treatment of some of the line excitation processes 
(Baldwin et al. 2004), resulted in a better understanding of 
several strong Fell features but changed little with regard to 
the main conclusion: there are several mechanisms capable 
of enhancing Fell emission relative to the hydrogen Balmer 
lines and increased iron abundance is only one of them . 

The correlation shown in Fig. [5] provides a new way to 
solve the iron abundance problem. As shown in S04, the N/C 
abundance (and metalicity in general) as measured from the 
NvA1240/CivA1549 line ratio (Hamann et al. 2002 and 
references therein) is strongly correlated with L/LEdd over a 
large luminosity range. This result is in conflict with the ear- 



lier studies by Hamann and Ferland (1993; 1999) who sug- 
gested that the N/C ratio depends on source luminosity (see 
also Dietrich et al. 2002; Dietrich et al. 2003 for the possible 
dependence on M(BH)). The S04 finding is based mostly on 
the high metalicity of low luminosity narrow line Seyfert 1 
galaxies (NLSls) that do not fall on the original Hamann and 
Ferland (1993) correlation (Shemmer and Netzer 2002). The 
high metalicity of NLS 1 s gains further support by the detailed 
analysis of the absorption line spectrum of Mrk 1044 (Fields 
et al. 2005). 

Given the dependence of metalicity on L/LEdd, and the cor- 
relation of FeII/H/3 with L/LEdd shown in Fig. [5] we suggest 
that FeII/H/3 is an iron abundance indicator for all sources 
with L/LEdd ~ 0.03. This conclusion does not rely on uncer- 
tain atomic models and complicated radiative transfer in pho- 
toionized BLR clouds. It is based on observed correlations 
in large AGN samples. There is no accurate way to calibrat- 
ing the iron-to-hydrogen abundance because of the uncertain 
modeling. However, we can compare this relationship to the 
S04 correlation of N V X 1240/C IV A 1549 vs. L/L Edd (their 
Fig. 6). The comparison shows that the slopes of the two 
correlations are similar. For example, the use of the BCES 
method with the assumption of a uniform error of 0.1 dex 
in both axes gives Fell/Hfi°cL/LEdd°' 7 ■ Assuming that Fe/H 
scales with O/H, and using the Hamann et al. (2002) calcu- 
lation to calibrate N/C vs. metalicity, we find that solar Fe/H 
corresponds to FeII/Hj3~ 0.6. 

3.3. The growth time of massive black holes 

Next we assume that BH growth is dominated by gas ac- 
cretion and BH mergers are of secondary importance. Given 
this assumption, we can calculate the growth time using the 
Salpeter (1964) original expression as 

fgrow . 4xl0«^l log -^-^yr, (6) 

L l L Edd M see d /active 

where / ac tive is the fraction of time the BH is active. For sim- 
plicity we assume a mass independent accretion efficiency of 
tj = 0.2, to reflect the assumption of a non-zero angular veloc- 
ity BHs, and a relatively large seed BH mass of M seed = 10 4 
Mq. Obviously there is a range in both properties (e.g. Heger 
and Woosley, 2002; Ohkubo et al. 2006) but the analysis is not 
very sensitive to either of those provided M seed is not much 
larger than assumed here. 

Given the above growth time, and assuming, /active = 1, we 
can calculate for each BH the time required to grow to its ob- 
served mass relative to the age of the universe at the observed 
redshift, t gr0 w/tuniverse- This is shown in Fig.[6]for all sources 
as functions of both L/LEdd an d redshift. Evidently, most 
(66%) BHs did not have enough time to grow to their mea- 
sured mass given their observed accretion rate. The fraction of 
such sources found here is, of course, an upper limit, since all 
BHs missing from our sample, due to their very small L/LEdd 
and the SDSS flux limit, also belong in this group of objects. 
Given the assumption of growth by accretion, we find that in 
the majority of cases, there must have been one or more past 
episodes where the accretion rate was higher than the value 
measured here. As shown in §4 below, this conclusion is of 
major importance in understanding the time dependent metal 
enrichment in AGNs. 

3.4. Redefining narrow line type-I AGNs (NLAGNls) 

The subgroup of NLSls is usually assumed to contain 
AGNs with the largest accretion rates (e.g. Grupe et al. 2004 
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FIG. 6. — BH growth-time, relative to the age of the universe, as a func- 
tion of L/LEdd (left) and redshift (right). Objects above the horizontal line at 
t(grow)/t(universe)=l must have had past episodes of faster accretion. Sym- 
bols are the same as those of Fig. 1 . 



FIG. 7. — A new defmion of narrow line type-I RQ-AGN (NLAGNls) 
based on accretion rate. The large symbols represent very luminous AGNs 
from S04 with a median BH mass of 10' 5 M®, The dividing line is drawn at 
L/L E dd=0.25. Symbols are the same as those of Fig. 1. 



and references therein). Such sources were first discussed in 
Osterbrock and Pogge (1985) who introduced the defining cri- 
teria of FWHM(HJ3)< 2000 km s -1 . NLSls are known to 
have other extreme properties such as extremely steep soft X- 
ray continuum and larger than average FeII/Hj3 . 

The recent study by S04 includes H and K-band spec- 
troscopy of about 30 high redshift, very high luminos- 
ity AGNs. In most of these sources L/Le c / c /> 0.4 yet 
FWHM(H/3) is much larger than 2000 kms" 1 . These au- 
thors suggested that the defining criteria for NLSls should be 
mass and luminosity dependent thus reflecting the fundamen- 
tal physics affecting the line width. This approach is different 
from the one taken by McLure and Dunlop (2004) who as- 
sumed the boundary between the groups of narrow and broad 
line type-I AGNs is at 2000 km s ~ 1 , independent of luminos- 
ity. Given this assumption, they did not find a larger fraction 
of NLSls among the highest luminosity sources in their sam- 
ple. 

Our SDSS sample contains a large enough range in both 
M(BH) and L/LEdd to study the S04 suggestion in more de- 
tail. To illustrate this we plot in Fig. accretion rates as a 
function of FWHM(Hj3) for all sources and shaw, on the same 
diagram, the two mass groups introduced earlier and 3 1 high 
redshift sources from S04. The median mass for the third 
group is about 3 x 10 9 Mq and the median L5100 about 10 46 - 6 
erg s _1 (an orders of magnitude larger than the most lumi- 
nous z < 0.75 SDSS sources). The diagram shows also a line 
at L/LEdd=0-25 that intersects the various subgroups at dif- 
ferent line widths. At small z, this line is consistent with the 
Osterbrock & Pogge (1985) definition of NLSls in terms of 
continuum luminosity and BH mass. We suggest that this line 
is a more physical measure of a general group of "narrow line 
type-I AGN" (hereafter NLAGNls). Using eqns.^and 2 with 
/l = 7, and normalizing the relationship to have a boundary 
at L/Le^=0.25, we get either a luminosity -based definition, 



FWHM(H beta ) < 1870 



^5100 



10 44 ergs 



-1 



0.175 



kms 



or a definition based on black hole mass, 



FWHM(H beta ) < 1700 



M 



10 7 M,. 



0.175 



kms 



-1 



(7) 



(8) 



3.5. FWHM([0 ill] A5007) as a black hole mass indicator 

The discovery of a tight correlation between black hole 
mass and the stellar bulge velocity dispersion (a*) resulted 
in several attempts to investigate secondary mass indica- 
tors related to emission line widths (cA This is particu- 
larly important in those cases where the strong non-stellar 
continuum prevents the measurement of a*. For example, 
FWHM([0 III] A5007) has been suggested as a proxy for cr* 
in type-I AGNs (e.g. Nelson 2000; Shields et al. 2003; Grupe 
& Mathur 2004; Boroson 2003; Boroson 2005). 

Boroson (2003) used a small subsample of 107 ob- 
jects with z < 0.4 from the early release SDSS data to 
study the FWHM([0 III] A5007) vs. ff* relationship. The 
FWHM([Om] A5007) in this study was measured from the 
observed line profile with no profile fitting. He finds a real 
correlation between FWHM([Olll] A5007) and the BH mass 
with a flattening of the relationship at higher luminosity and 
a very large scatter (a factor of ~ 5). In a later publication, 
Boroson (2005) revisited the same problem using a larger 
sample (about 400 sources) covering the same redshift inter- 
val and focusing on the systemic velocity indicated by the 
[O III] X50Q1 line center. He finds a systematic blueshift of 
the line peak with a magnitude which is correlated with the 
line width and with L/LEdd- An extensive discussion of the 
merit of such techniques, including an analysis of a large sam- 
ple (1749 object) of type-II AGNs, with a median redshift of 
0.1, where both a* and the width of several narrow emission 
lines have been measured, is given in Greene & Ho (2005). 
According to these authors, FWHM([0 in] A5007)/2.35 is a 
good proxy for <7» provided the strong line core, without the 
commonly observed extended blue wing, is used to define 
FWHM([0 III] A5007). Greene & Ho further suggested that 
deviations between a* and FWHM([0 III] A5007)/2.35 are re- 
lated to L/LEdd- Large deviation probably reflect non-virial 
motions due to a strong radiation field. The [O III] A 5007 line 
luminosity has also been used as a primary bolometric lumi- 
nosity indicator for type-II AGNs (e.g. Heckman et al. 2004) 
whose non-stellar continuum cannot be directly observed. As 
shown by Netzer et al. (2006), this is associated with a large 
scatter, and perhaps also a luminosity and/or redshift depen- 
dence. 

Our SDSS sample is different from all earlier ones by being 
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FIG. 8.— The correlation of M(BH) vs. FWHM([0 III] A5007) for two 
redshift bins: z=0. 1-0.2 (bottom) and z=0.4-0.6 (top). The entire RQ sample 
is plotted as small points in both panels and the objects belonging to the 
specific redshift group are plotted as large black symbols. The two lines in 
each panel are best fit curves for the redshift group in question (solid line) 
and the Tremaine et al. (2002) fit to the M(BH)-o~, relationship (dashed line). 
The two fits are basically identical for the z=0. 1-0.2 group. 



much larger and by covering a wider redshift range. It can 
thus be used to test the earlier findings regarding type-I AGNs 
at higher luminosities and redshifts and with a better statis- 
tics. Unlike Boroson (2003; 2005) and Greene & Ho (2005), 
we do not restrict ourselves to high S/N spectra and attempt 
to fit all sources with a detectable [O III] A5007 line. As ex- 
plained in §2.2, our FWHM([0 III] A5007) is based on a sin- 
gle Gaussian fit to the core of the line, similar to the method 
of Greene & Ho (2005). We have correlated these measure- 
ments with M(BH) and L/LEdd measured here and show the 
results in Fig||]as M(BH) vs. FWHM([Olll] A5007) for the 
entire sample as well as for two redshift groups: z=0. 1-0.2 
and z=0.4-0.6. To guide the eye, we have calculated the best 
slope for each of the groups using the BCES method assum- 
ing a uniform error of 0. 1 dex and plot them on the diagram 
together with a curve representing the Tremaine et al. (2002) 
parametrization of the M(BH)-(7h. relationship. 

Fig. [8] clearly shows the non-linear dependence of 
log(M(BH)) on log(FWHM([Offl] A5007)). Expressing the 
correlation as M(BH)°cFWHM([Olll] A5007)"W we find a 
range of a(z) as listed in Table 2. As evident from the table, 
the slope is decreasing with redshift and the correlation be- 
comes insignificant at z > 0.5. As a curiosity, the agreement 
with the Tremaine et al. (2002) best slope (a — 4.02) is al- 
most perfect over the redshift interval of 0.1-0.2 (see bottom 
panel of diagram). We suspect that the overall agreement be- 
tween a t and FWHM([Om] A5007)/2.35 found by Grupe & 
Mathur (2004) and by Greene & HO (2005), is the result of 
those samples being dominated by low redshift sources. 

High accretion rate can also influence the [Olll]A5007 



TABLE 2 

Fit coefficients for M(BH) vs. FWHM([Oiii] A5007) as defined 





IN 


§3.5 




Az 


«(z) 


L/ L Edd 




0.05-0.15 


5.8± 1.5 


0.01-0.03 


3.4±0.1 


0.10-0.20 4.0 ±0.6 


0.10-0.15 


2.7 ±0.7 


0.37-0.43 


1.6±0.2 


0.30-0.50 


2.1±0.2 


0.66-0.75 


1.5±0.9 


all 


2.9 ±0.4 



line width by introducing a non-virial component to the 
line profile, perhaps due to acceleration by radiation pres- 
sure force. This has been studied by Netzer et al. (2004), 
Boroson (2005) and Greene & Ho (2005). Netzer et 
al. found that FWHM([Olll] A5007) correlates with L5100, 
M(BH) and L/LEdd but the dependence on L/LEdd is weaker 
than the dependences on the luminosity and the mass. To 
check this in the present case, we divided the SDSS sam- 
ple into groups according to L/LEdd and assumed that 
M(BH)ocFWHM([0 III] X5QQl f {L l L ^. We find a clear de- 
crease of /3 with L/LEdd- The computed values are given 
in Table 2. Thus, larger L/LEdd results in a flatter M(BH)- 
FWHM([0 III] A5007) relationship. 

The redshift dependence of the M(BH)- 
FWHM([0 III] A5007) relationship can be explained in 
several ways. First, selection effects must be important 
as more and more low mass, low luminosity and small 
L/LEdd sources drop from the sample at higher redshifts. 
We suspect that the paucity of small mass black holes at 
high redshifts cannot explain the entire change as the overall 
range in M(BH) at z > 0.5 is quite large (about 1.3 dex). 
The increasing number of high L/LEdd sources is probably 
more important as this tends to flatten the relationship. This 
is in line with the Greene & Ho (2005) suggestion and was 
tested quantitatively by using their recipe (their eqn. 1) to 
correct the observed FWHM([0 III] A5007) for L/L E dd- The 
procedure results in a somewhat better agreement with the 
Tremaine et al. (2002) relationship over part (but not all) 
of the range in FWHM([Olll] A5007). This by itself does 
not explain the observed spread. We conclude that on top 
of a large intrinsic spread in all quantities related to the 
[O III] A 5007 line, there are also systematic effects that limit, 
severely, the usefulness of this approach. Thus, the methods 
used to obtain bolometric luminosity and BH mass from the 
luminosity and the width of the [O III] A 5007 line suffer from 
a large scatter and, perhaps, systematic uncertainties. 

An interesting idea is a real redshift evolu- 
tion of the M(BH)-(Xh relationship. In this case, 
FWHM([Olll] A5007)/2.35 can be used as a surrogate 
for a* at small z yet the non-parallel evolution of the bulge 
and the BH mass requires a modification at larger redshifts. 
Further discussion of this idea is beyond the scope of the 
paper. 

3.6. Other correlations 

The SDSS sample is an invaluable source for studying sev- 
eral other correlations of line and continuum properties in 
type-I AGNs. Such a study goes beyond the scope of the 
present work that focuses on mass, accretion rate and met- 
alicity as a function of redshift. Since most those correlations 
can easily be computed, we list them below, with only a brief 
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TABLE 3 

Correlations with a chance probability smaller than 10~ 5 



Property A 



Property B RQ-AGNs RL-AGNs M(BH)=10 7 - 5 ' 1 M, 
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FeII/Hj3 


L5100 
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FeII/Hj3 
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L/L Edd 
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FeII/Hj3 


EW(Hj6) 








Hj3(red)/H/3 (total) L 510o 
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Hj3(red)/Hj3(total) M(BH) 
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+ 




Hj3(red)/H/3 (total) EW(H/3) 


+ 


+ 


+ 


Hj3(red)/Hj3 (total) FeII/Hj8 









discussion and leave the more detailed investigation to a forth- 
coming publication. 

Table 3 shows the results of our Spearman rank correla- 
tion analysis for several line and continuum properties. Sig- 
nificant correlations, with a chance probability smaller than 
10~ 5 , are marked by either a + or a — , depending on the slope 
of the correlation. In particular, we find a negative Baldwin 
relationship (Baldwin 1977) for both Hj3 and the Fell lines. 
The former has already been found by Croom et al. (2004) 
in their study of the 2dF AGN sample. We also find a ten- 
dency for the fractional luminosity of the red part of the H/3 
line (measured relative to the position of [O III] A5007) to in- 
crease with continuum luminosity. This is an opposite trend 
to the well known blue shift of the C IV A 1549 line in high 
luminosity sources and does not seem to be associated with 
L/LecUI- We note that the definition of the systemic velocity is 
a bit ambiguous since we rely on the measured wavelength of 
the [O III] A5007 line. Yet, this lines was shown by Boroson 
(2005) to have a noticeable blueshift in large L/LEdd AGNs. 

4. DISCUSSION 

The main new results presented here are the measurement 
of the accretion rate, in several mass groups, as a function of 
redshift and the correlation between black hole growth rate 
and BLR metalicity. 

4. 1 . Black hole growth and BLR metalicity 

The new results concerning the iron abundance and the 
growth time can be combined to infer the changes in the 
metal abundance of the BLR gas as a function of time. Since 
FeII/Hj3 is an iron abundance indicator which depends on 
L/LEdd, an d since most type-I AGNs had higher accretion 
rate episodes in their past, we cannot avoid the conclusion that 
those past episodes were characterized by enhanced metalic- 
ity relative to the one measured here. Computing BH growth 
times for the 3 1 high luminosity sources in S04, and noting the 
dependence of N/C on L/LEdd found by those authors, results 
in a similar conclusion, despite of the much larger redshifts 
of these sources. Thus, BLR metalicity can go through cycles 
and is not monotonously decreasing with decreasing redshift. 

The emerging picture of a complex metalicity, accretion 



rate and cosmic time dependences is not easy to explain. One 
possibility is that episodes of fast accretion onto the central 
BH are associated with enhanced, large scale star formation 
episodes that lead to increased metalicity throughout the host 
galaxy or at least in its central kpc. If the enriched mate- 
rial finds its way to the center of the system, increasing both 
L/LEdd an d the metal content of the BLR, it would explain 
all observed properties. Such a scenario must also involve 
ejection and/or accretion of the enriched gas toward the end 
of the active phase such that the metal content of the BLR at 
later times, when the fast accretion phase is over, is reduced 
relative to the peak activity phase. The formation of enriched 
BLR clouds from the wind of an enriched accretion disk is 
just one of several such possibilities. 

The (somewhat uncertain) dependence of FeII/H/3 on red- 
shift suggests that the general, time averaged BLR iron abun- 
dance is slowly increasing with cosmic time while undergo- 
ing much larger fluctuations associated with the fast accretion 
rate episodes. Finally, if the time scale of channeling such en- 
riched starburst gas into the center is short compared with one 
episode of BH growth, we would expect to see strong starburst 
activity associated with the BH activity in many AGNs. 

4.2. Black hole growth and cosmic star formation 

Our new results clearly show the changes in BH mass and 
accretion rate with redshift. For example, the histograms 
shown in Fig.[5]clearly illustrate that the smaller BHs are the 
fastest accretors at all redshifts studied here although large 
mass BHs at high z can accrete as fast as small mass BHs at 
low z. 

Given the results shown in Fig. |4] and listed in Table 1, we 
can compare the BH accretion rates found here and the known 
star formation rate over the redshift interval of 0-0.75. For 
this we use eqns.|4]and[5]and various publications describing 
evolutionary models and star formation rate measurements 
(e.g. Chary & Elbaz 2001 and references therein). A simple 
fit to some of those results gives p t ^ £/ or p» « (1 +z) s where 
1.4 < 7 < 2 and 6 < 8 < 9. The large range of slopes repre- 
sents the uncertainties on the data as well as on the various 
evolutionary models (see Chary & Elbaz 2001 for more infor- 
mation). This is a universal rate averaged over all starburst 
galaxies of all morphologies and luminosities. The range of 
slopes is somewhat narrower than the one found here for BH 
growth rates (Table 1) and the middle of the range is in rough 
agreement with the growth rate of M=10 8 M© BHs. Thus, the 
average BH growth rate for z < 0.75 AGNs seems to agree 
with the star formation rate over the same redshift interval. 

Marconi et al. (2004) used a detailed modeling of sev- 
eral observed AGN luminosity functions to argue that ac- 
cretion onto BHs, integrated over the entire AGN popula- 
tion, proceeds at a rate similar to the star formation rate at 
small redshifts. Their model fitting to the data over the red- 
shift range 0-0.6, assuming L/Lsdd = ^< can be described by 



L(all AGN) 



(1 



\4.2 



(A. Marconi, private communica- 



tion). Obviously, the Marconi et al. (2004) results cannot 
be directly compared with ours. Their total accretion rate (i.e. 
luminosity) was calculated for the entire AGN population and 
the redshift dependence was obtained under the assumption 
of a uniform L/LEdd- 

The new results shown here go one step further by neglect- 
ing the uniform L/LEdd assumption and by demonstrating that 
different mass groups are characterized by different growth 
rates. A more detailed comparison between the two stud- 
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ies requires an additional information about sources that are 
missing from the SDSS data but are found in X-ray selected 
samples. This is beyond the scope of the present work. 
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